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Abstract 

The exact large deviation function (Idf ) for the fluctuations of the energy density field is computed 
for a chain of Ising (or more generally Potts) spins driven by a zero-temperature (dissipative) Glauber 
dynamics and sustained in a non trivial stationary regime by an arbitrary energy injection mechanism 
at the boundary of the system. It is found that this Idf is independent of the dynamical details of 
the energy injection, and that the energy fluctuations, unlike conservative systems in a nonequilibrium 
state, are not spatially correlated in the stationary regime. 
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1 Introduction 



In a series of recent papers, Derrida et al. ^ studied the noncquilibrium stationary states of exactly 
solvable models characterized by conservative inner dynamics: particles diffuse in a one dimensional 
chain of sites without annihilation, and are described at the hydrodynamical length- and time-scales 
by a Fick's (diffusion) law (with or without a systematic drift force). The authors succeeded in exactly 
computing the noncquilibrium free energy (associated with a given density profile) which displayed very 
interesting features: for instance, the noncquilibrium situation induces long range correlations which 
make the free energy non additive. The authors discovered however an elegant and quite unusual 
sub-additivity principle whose generality for other models is still an open question. Their results also 
showed that although the second moment of density fluctuations is the same as previously found with 
fluctuating hydrodynamics arguments a quadratic (gaussian) approximation is not sufficient for 
higher cumulants. 

In this paper, we address the same question of characterizing of noncquilibrium stationary states 
(NESS), but for dissipative inner dynamics. Such situations are widely encountered: for instance, the 
granular matter continuously shakcd by a piston can be considered for long time scales as being in a 
noncquilibrium stationary state; the velocity field of a turbulent flow is also subjected to a dissipative 
(Navier-Stokes) equation and can be sustained in a stationary state as well. These dissipative systems 
are altogether different from conservative ones, since no concept of thermodynamical equilibrium can 
be applied to them: in absence of energy input, the full rest is the ultimate state of the system; 
as a result, no pertubative strategy can be deployed to describe their physics, and in fact the non 
equilibrium stationary states of dissipative systems are not well known (see for instance 0] and 
references therein). 

In particular, at the light of the recent works aforementioned, natural questions arise for dissipative 
NESS: are they similarly characterized by a high level of correlation ? are the fluctuations to some 
extent insensitive to the details of the inner dynamics and/or the injection mechanism ? is the notion 
of noncquilibrium "thermodynamical potential" even relevant ? 

To answer these questions, and to develop further the ideas presented in we consider here a one- 
dimensional dissipative model for which some exact results on the energy structuration can be extracted: 
our model is a semi-infinite ID chain of Ising spins (easily generalized to a Potts model) subjected to 
zero-temperature Glauber dynamics. These intrinsically dissipative dynamics are supplemented with 
an arbitrary (Poissonian to be simple) flipping process of the first spin as a way to inject energy into 
the system. 

Our results can be summarized as follows: i) the concept of noncquilibrium potential can be extended 
via a non homogeneous space coarsening of this dissipative system; ii) this potential is independent of 
the nature of the injection mechanism, and reflects mainly the dynamical inner self organization of the 
system; iii) this potential does not display correlations, i.e. it is "additive" , as soon as a stretching of 
the space is taken into account; iv) such a potential cannot be defined for any dissipative system, and 
the conditions that system dynamics must fulfill to have a global observable obeying the large deviation 
theorem are not clear. 

The paper is organized as follows. In a first part we define the model and its dynamics; then we 
give some preliminary results concerning the stationary state and introduce the cumulative energy E. 
Afterwards we compute the large deviation function of E and the Idf associated with an energy profile. 
We end up with a discussion of the results. 

2 Definition of the model 

The model we consider is a I-dimensional semi- infinite chain of Ising spins ao, ai, . . . , aj , . . . (a finite 
chain could be considered as well). The spins follow zero-temperature Glauber dynamics, except for 
the spin cfq , which flips according to a Poisson process of parameter A : the probability associated with 
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a flip of spin j from a conflguration C between t and t + dt is w{C — > Cj)dt = [1 — ^^(ctj+i +(Tj_i)/2] x 
if j ^ (the starting configuration is termed C and the j-flipped configuration Cj) and Xdt for tro- As 
a result, the master equation describing dynamics in the model is 

dtP{C) = - ^ w{C ^ C,)P{C) + J2 ^ C)P{C,) (1) 

= -A[P(C) - P(Co)] - - (7,(a,+i + a,_i)/2]F(C) + ^^[l + <T,(a,+i + (7,_i)/2]P(C,) 

(2) 

In the last equality, the spin values refer obviously to configuration C. 

Inner Glauber dynamics are essentially dissipative: the domain walls (which are elementary energy 
excitations) move randomly and annihilate by pairs when colliding. Thus, in absence of energy input, 
any initial condition would eventually relax to a state characterized by the same value of all spins; in 
the model considered here, such an external input is provided by the Poissonian motion of the spin 
labeled 0, which gives energy to the system when ctq fiips from tJi to —ai. As a result, after a transient, 
a stationary state takes place which is described in the following^ . 

3 The stationary state 

3.1 Mean injected power 

If the summation croCi(- • ■) is made on the dynamical equation one gets 

dt{croo-t) = -2(A + l){aoai) + (o-qo-j+i) + (cro(T.i_i) (3) 

When the stationary state is assumed, and according to o-g = 1, (O leads to (o'oo'j) = [1 + A — \/A^ + 2A]' 
(this correlation vanishes when i — s- oo). 

From this calculation an interesting physical result can be deduced: the (mean) power Pinj injected 
by the spin inside the system is related to (aoai), since between t and t + dt, the energy ceded is in 
average 2Xdt x [Prob(cro = ci) — Prob(cro = — ci)], whence [7] 

Pi„j = 2A(CToai) (4) 
= 2A(1 + A - VA2 + 2A) (5) 

The injected power is an increasing function of A which saturates to 2 for large A; this value is phys- 
ically dictated by internal dynamics, i.e. the ability of the system to diffuse into the bulk the energy 
excitations created at the boundary. It is to note that there is thus no notion of "optimal" time scale 
concerning the energy injection, which is not evident a priori: one could imagine that an optimal 
waiting time could leave domain walls move away from the boundary. 

3.2 Two-points correlations and energy density profile 

Another interesting quantity describing the stationary state is the average energy density profile 

(e„) = 1 - (cr„(T„+i) (6) 

This quantity is more complicated to obtain. To this end, it is useful to interpret Glauber dynamics 
as coalescing paths dynamics |S]: the update of any spin is equivalent to a random choice (Poissonian 

^we are aware that the stationary state of the infinite system needs in general a diverging time to be established, but this 
problem can be circumvented by a clever choice of the initial conditions 
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with parameter 1) of the spin value among one of its two neighbours, such that the value cr„(i) can 
be traced back in time from spin to spin until either the time origin or the boundary ctq is reached 
at a certain spin index or time (respectively). If the stationary state is studied, the notion of time 
origin is irrelevant, and any "path of constant spin value" eventually reaches the zeroth spin (see fig. 
If two spins (T„ and Up are now monitored, two paths emerge from n and p, are traced back in 



-t 




j spin index 

Figure 1: The value of (Tj{t) can be traced back in time from neighbouring spin to neighbouring spin until 
it reaches ctq whose statistical properties are known. 

time, and possibly coalesce if they meet each other before reaching the boundary. They are statistically 
independent, since the flipping processes in two distinct sites are not correlated. As a result, the average 
{(JnO'p) can be calculated as follows: consider two random walkers starting at sites n and p. If they 
meet before touching the site 0, this occurence of the paths gives a factor 1 in the computation of 
(CT„(Tp). Otherwise, the n-path (resp. p) arrives at site zero at time i„ (resp. tp), and this occurence 
gives a factor (o'o(^n)'''o(^p)) = exp(— 2A|i„ — tp|), where the average is on the Poissonian process ctq. 
As a result original dynamics are mapped on dynamics of two random walkers plus dynamics of ctq. 
These considerations can be summarized as 

/•oo pOC 

+ / dtn / dtpProh[tn,tp and {n,p)]cxp{~2X\tn - tp\) (7) 
Jo Jo 

In this expression, c„p is the probability that the walkers starting at n and p meet each other before 
reaching the zeroth site (we will also use Cn,p ~ 1 — c„p; generally, a comma is put between two 
indices when the associated walkers are supposed to avoid each other; conversely the absence of comma 
holds for coalescing walkers); besides Prob[i„,tp and {n,p)] is the probability that the n- walker and 
the p-walker reach site at times tn and tp respectively without having met each other beforehand. 

These expressions can be exactly computed. Indeed, consider two independent random walkers 
starting at n and p. There are a priori three possible situations (fig. 13): either (i) they do not cross 
at all, or (ii) they cross and "exchange" their arrival time, or (iii) they cross without exchanging their 
arrival time. It is easy to see that each situation of type (ii) can be mapped to a situation of the type 
(iii), and vice- versa (fig. Moreover, the probability is of course conserved by the mapping, since 
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(i) (ii) (iii) 

Figure 2: Three possible generic situations for two independent walkers 



the two paths are of the same statistical nature. As a result, if n < p: 



Cn,p = Prob(path) — Prob(path) 

type (i) and (iii) paths type (ii) paths 

din I dtp [Pn(tn^Pp(tp^ -Pnitpj-Ppitn}] 







^np ~ 2 



{t-p')Pp (tji ) 



tr^<t„ 



(8) 

(9) 
(10) 



where Pn{T) is the probabiUty (density) for a walker initially {t = 0) at n to reach the site for the 
first time at i = r. We get, for n ^ p: 

{<7n<Jp)^2 dt dt' Pn{t')Pp{t) + dt dt' [Pn{t)Pp{t')- Pn{t')Pp{t)]cxp{~2X\t~t'\) 

Jo Jt Jo Jt 

(11) 

We are going to see that the dominant contribution to {(T„(7p) in the limit ?i ^ oo is given by the first 
term only. This fact has important consequences on the structure of the stationary regime. 



The probability P„(t) is obtained as a solution of a Brownian motion with an absorbing boundary 
at site 0: the probability p(m,t\n) that a Brownian walker starting at site n at t — he at site m at 
time t without having touched the zeroth site obeys the dynamical equation dtp{m\n) = p{m + l|n) + 
p{m — l|n) — 2p(m|n) if m > 1 and dtp{l\n) = p{2\n) — 2p{l\n). The probability density P„(t) is simply 
given by Pn{t) = p{l,t\n) (actually Pn{t) is p{l,t\n) times the transition rate associated to the jump 
1 — > 0, but this rate were chosen equal to one). That p{m, t\n) = Pf{m — n,t) — Pfijn + n, t) (with pf 
is the free-boundary Brownian motion starting at site 0) is the solution is directly verified, whence one 
gets 

Pn{t) = 2 e^^* sin(27rg)sin(27r(jn) (12) 

Jo 
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Figure 3: mapping of type (ii) to type (iii) paths 



With this expression, one can recast c„p in 

,1/2 



4 y dq cot(7r(7) sin(27rqri) ^2 sin^ nq + l — 2\j wv" ixq + sin^ Trg 



An asymptotic expansion leads to the foUowing resuh, valid for n ^ cx),p > n: 



-At£ 



The second term of can be estimated using the asymptotic behaviour of P, 
1 



Pninr) 



1 



n-,oo^27rn r(4T2 + 1)1/4 



exp n 



1 - 2t - log 



1 
27 



1 

4^2 



1 



r»iy4^ 



:exp(-^) 



A cumbersome computation yields for n oo and p > n: 



dt 



dt' [Pn{t)Pp{t') ~ Pn{t')Pp{t)] exp(-2A|i - t'\) 



24 ri,p(p^ 



7rA2 (n2+p2)4 



(13) 
(14) 

(15) 
(16) 

(17) 



It is important to remark that in the asymptotic limit n —f oo, the dominant contribution to ((JnO"p) 
is given by the "coalescing term" c„p, which is independent of the value of A ; moreover the dominant 
contribution is (besides the constant 1) non summable and leads to a logarithmic divergence, whereas 
the A-dependent corrections are of order (for n — p ~ 0(1)) and therefore obviously summable. 
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These results can be applied to the energy density profile: 

(e„>c^ — + of4) (18) 
Once again the leading term of the energy density profile is independent of A and non summable. 

3.3 Energy density correlations 

Once we have the averaged profile (e„), it is natural to subsequently ask for the correlation function 
of e„: 

The 4-point correlator (cr„(T„+icrpCrp+i) can be written (cr„(T„+icrpCrp+i) = c„„+ipp+i + c„„+i,pp+i + 
x{n,p) where c„„+ipp+i is the probability that four walkers starting at n, n + 1, p and p + 1 coalesce 
before reaching the boundary; c„„+i^pp+i (note the comma) the probability that the walkers "n" and 
"n + 1" on one side, walkers "p" and "p + 1" on the other side coalesce and that the two groups reach 
the boundary separatly; and xi'^^iP) is the term arising from situations where one or several walkers 
reach the boundary without coalescence. This x("-5P) is obviously A-dependent and is composed of 
several terms, each corresponding to a particular scenario for the walkers: with natural notations, 
{n,n+l,p,p+l}, {{n n+1 p),p+l}, {n,{n + l p p+1)}, {{n n + l),p,p+l}, {n,{n + l p),p+l}, 
{n, (p p+1)} are the different possible situations where A plays a role. Let us consider the first 
term: 



{n,n+ l,p,p+ 1} = / Prob[ii,t2,t3,t4 and (n, n + l,p,p+ l)]((To(ti)fTo(t2)o-o(t3)o-o(i4)) 

Jti<t2<t3<t4 

(20) 

(the probability term refers to four walkers starting from the four sites considered and reaching without 
coalescence the site at times ti. . .t^ respectively). This integral must be approximately evaluated; 
its dominant term comes from regions in the (^1,^2,^3,^4) plane where ti ~ t2 and ts t^: in these 
regions only is the 4-point correlator in ctq significantly different from zero. As a result, 

{n,n+l,p,p+l} - A"2 / Prob[ii,ti + A"\i2,i2 + A~^ and {n,n + l,p,p + 1)] (21) 

Jti<t2 

<X-^ / Prob[ti,ti + A"^ and (n,n + 1)] X A"^ / Prob[t2, ^2 + A"^ and(p,p+l)] 

Jti Jt2 

(22) 



But as 

A-i 



d< Prob[t,< + A-^ and (n,n + l)] - A'^y dtP„(t)Pn+i{t + X''^) - Pn(t + X-'^)Pn+i{t) (23) 



we conclude that {n, n + 1, p, p+1} gives a very rapidly decreasing contribution, and we can notice 
besides that X]n<p{"'' "^^ -i'P'?'+ -i} ^ "^^^ situation is similar if other A-dependent terms are 
considered. For instance, 

{(n n-fl p),p+l} - A^W dtProh[t,t + X^^ and {{n n + 1 p),p+l)] (25) 

"'0 

< X-^ / dtProh\t,t + X-^ and (n,p+l)] - X'^pn-^ (26) 

Jo [P^ + n^r 
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(28) 



(note the similarity between this approximation and the exact Umit of H17|l ). Once again this term is 
very small and verifies X]ri<p{("' "- + 1 p)iP+^} < oo- 

A special attention must be paid to the term {n,n + l, {p p+1)}. Actually it is not summable: 
n + 1, {p diverges with a contribution 0{N) and another 0(log A^). But this term 

is partially "disconnected" and these divergences are exactly balanced in (|19|l by a term coming from 
(cr„cr„+i) (cTpCTp+i) = (c„„+i + {n,n+l})(cpp^i + {p,p+l}): actually {n,n+l}cpp^i is the counterterm 
which kills these "hybrid" divergences. 

So, the correlation function of the energy can be recast as 

1 + C„,„+l + - Cn,n+lCp^p+l + x("-;P) (27) 

where x(^iP) depends on A and verifies X]n<p x("iP) < We will see in the following that 

Cnn+lppH-l "1" Cnn+l,pp4-l 1 ^n.p}-l + ^n.p + ^rti-l.p}-! ^n.n+l ^rt^l.p ^p^p}-! 

~t~ f^n,n+l^p.p^l + Cn.p{-lCnr{-l.p ^p^rh\-l ,p\-l 

whence finally 

{^n^p} ('^n) ^n,p^l + ^n,p + ^n+ljjH-l ^n+l.p + ^7i.pj-l^n+l,p ,p^rt+l .p{-l + X(^i7^) (^9) 

Local correlations are then easily deduced from (|14|) : 

32n(n+l)(2n+l) 1 
(e„e„+fc) - (e„)(e„+fc) ~ — (30) 

n largo .fe^oo OTT'' K'^ 

I \ I M \ 32(n+l)(2n+l) 1 

n large ,k^oo AtT Tl^ k° 

and (e^) - (e„)2 ~ 4/(7m). Thus, the energy distribution has anticorrelations which decrease with the 
distance. The second scaling is clearly much more relevant to describe the system, a fact that will 
become clearer in the following. 

3.4 The cumulative energy E 

It has been seen that the local structure of the stationary energy field is of course affected by the in- 
jection mechanism. Nevertheless the preceding preliminary computations highlighted a very important 
"boundary layer" phenomenon. To precise this, let us define the cumulative energy E: 

(32) 

ri=0 

which is nothing but the energy of the A^ first spins. The average of E and its fluctuations share the 
same following properties: 

• they are 0(log A) when N ^ oo 

• their leading term is independent of the injection mechanism. 

To sec this, let us remind first that (E) ^ ^ log A^ (from cq. (|18|l 'l. As to the fluctuations, the preceding 
section allows to write 

AE'^ = (E^) - {Ef (33) 



2 [—Cn.p+l+Cn,p+Cn+l,p+l — Cn+l^p + Cn^p+lCn+l^p — Cn^pCn+l,p+l\+ [(^n) ~ (^r^ 

(34) 



TV— >cx) 



l--)logAr (35) 
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The key point of this paper is to remark that these properties are actually shared by all cumulants of 
the distribution of E. We do not intend to give a demonstration for this, but owing to the arguments 
presented in the preceding subsection, this conjecture is reasonable: in an event of 2m walkers starting 
from sites located asymptotically far from the boundary, the probability associated with scenarios 
where one or several pairs of walkers reach the boundary is cither negligible or cancelled out by the 
cumulant structure of the mean. 



From the physical point of view, this property reflects a boundary layer structure: the injection 
mechanism provides energy into the system which spreads out quite efficiently despite the dissipation 
mechanism; however the dissipation is strongly efficient near the boundary and smoothes out the 
"memory" of the injection details. In other words the dissipation sets up an autosimilarity regime for 
the energy, which is independent of the amount transfered to the system. 

From the analytical point of view, this property has two important consequences. First, we shall 
see that there is a large deviation function (Idf) / associated with E, namely 

N^oo log TV 

and secondly this Idf is independent of the injection mechanism, for the cumulants of E arc directly 
related to the Taylor expansion coefficients of /. As already noticed [SHH, here again the Idf captures 
the essential features of an observable and fades away some "irrelevant" details: the summation process 
tends to "universalize" the different possible behaviours in some way. Moreover, this property allows 
the analytical computation of /(C), which is the subject of the subsequent section. 



4 The large deviation function of E 

As / does not depend on the injection mechanism, we will consider henceforth the limit A ^ oo 
(extremely quick flipping of (Tq). To get /, we will compute first the characteristic function G(/i) = 
{exp{— fj,E)) , which is the Laplace transform of the probability P{E); we will show that G(//) is dom- 
inated by a term exjp[g (fj,) log N]. Once g is known, the Laplace inversion formula shows that /(C) is 
given by a Legendre transformation: 

/(C) = max (^g{^i) + (37) 

provided that no analyticity breaking of the prcfactors of the exponentials in G occurs at the saddle 
point in the /i space (see |S] for instance). We will assume in the following that no such problem arises, 
which is the most often case. 

Using a classical trick of Ising spin systems and defining r = (exp(— 2/i) — l)/2, one can write G as 

G{fi) = ^ n e-A'd-.-.+i)^ = (j[ [1 + r(l - a.a,+i)]^ (38) 
= 1 + T^(l - (T,a,+i) + ^((1 - (7,a,+i)(l - ajaj+i)) + . . (39) 

i i<j 

= 1 + tSi+t^S2+t^S3 + ... (40) 

In the limit A — > c», all the averages in the preceding sum can be interpreted as probabilities related 
to walkers. For instance, (1 — CTidi+i) is just the probability c^.i+i that two walkers emerging from i 
and i + 1 do not meet before reaching the site 0; similarly ((1 — aiai^i){l — ajaj^i)) is (for i ^ j) 
the probability that during the wandering of four walkers emerging from sites labeled i, i + 1, j , j + 1, 
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walkers i and i + 1 from one side and walkers j and j + 1 from the other side have not collapsed (this 
interpretation is also valid if i + 1 = j). We will denote this probability pi^j in the following. We will 
see that the probabilities pi.j,K,... and therefore the terms Si can be expressed in terms of the Ci.j only. 



4.1 The first terms 

The term S'l is simply 5*1 = where the summation ranges from i = Otoi = A^ — 1. The 

term 5*2 can also be expressed in terms of the Cij. To this end, let us consider pij. The event (/, J), 
which is "walkers i and i + 1 do not collapse, as well as walkers j and j + 1" is equivalent to the event 
+ + 1) or {i, [i + Ij), j + 1", such that (remind that a comma holds for "have not met" 

whereas an absence of comma for "have met" ) 

Pi,j = Ci,i+i.],j+i + Q/H-iij+i (41) 

This probability can be transformed slightly. We have the important relations (actually gencralizablc 
to any situation) 

Ci,jk,e + Ci.j.M + Cij^kj = Ci^j^e ^^^^ 

^ij^ki ^ijj^k^i — ^i^j,k 

which are valid whatever i ^ j ^ k ^ £ (essentially they express the fact that if one enumerates 
correctly what can happen to an emerging path — all remaining paths being besides fixed — , and sum 
the corresponding probabilities, one gets the probability of the event whence that path is eliminated). 
They show how to express c... where some indices "stick" together, in terms of c... with no sticking. 
This yields for pj j and 6*2: 

Pi,j = Ci/j+i + Cis+ij+i — Ci^i+ij (43) 

'S'2 = ^ Ci^j+l j.j + l + ^ Ci^j+i,Ar (44) 
i<j i 

The complexity of the terms increases very fast with the number of indices: 

Pl,J,K = Ci,i+1 + + Ci^i+ijj+i^k,k+l + CiS+lJJ + lk,k+l + Ci,i+lj,j + lk,k+l (45) 

S3 = jj + i + Ci^i+i,jj + i^k,k+l + Ci,i+l,j,j + l,N (46) 

i<j i<j<k i<j 

(the sum S3 is obtained after a quite lengthy computation!). The next order requires some patience, 
and one gets 

Si = Ci,i+1 j"j'+l,fe,fe+l/,f+l + 2 Ci_i-|-i + Ci^i+ijj-i.i^k,k+l,N + Ci,i+1 j'j" + l,Af 

i<j<k<l i<j<k i<j<k i<j<k 

(47) 

No evident regularity emerges from these first terms. It is worth noting that a rewriting of these terms 
is possible, since it is possible to write a term c... with odd number of indices as a sum of c... with even 
number of indices, but this rewriting makes the situation all but more transparent. 

4.2 Resummation of G 

A resummation of G(/i) is nevertheless possible. To do this, it is convenient to introduce a schematic 
representation of the terms appearing in the Sj. With evident definitions, we introduce the notations 
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C(^') and C}j^ 

Si = C(i) (48) 
^2 = + dj^^ (49) 

^4 = Cf'^) + 2C(3) + Cj,^) + C^,^) (51) 

We will show that a recursion relation exists between S'„+27 Sn+i and Sn- The starting point is, besides 
the remark that Sn is a sum of terms C'^^' and C^^ , the following relation, obtained along similar lines 
as eq. (|^ : 

Pl,J,K =PlJ and j + l,fc,fc+l +PlJ and j + lfc,fc+l (52) 
= P/, J and j + l,A:+l ^ Pi, J and j + l,fc + P/, J and j + l,fe,fc+l (53) 

In this relation, pjj and j+i,k terms the probability that i and i + 1, j and j ' + 1, j + 1 and k do not 
meet pair by pair. The other terms have similar definitions. Moreover, we particularize the recursion 
relation on pi.j^K, but it is evident that a such relation is also true if some capital indices are implied 
at the left of /: this relation is true at any order actually. As a result, 

Pl,J,K ='^^Pl,J a.nd j+l,N + ^ P/.J and j + l,fc,fc+l (54) 
i<j<k i<j i<j<k 

Naively, we are tempted to deduce from this relation the formal law to pass from Sn to Sn+i- "put 
an extra TV on one side, increment the rank on an other, and add" . This rule is in fact not simple to 
implement for a recursion, since one does not know a priori how to interpret terms like (C^^)Ar and 
(C^''')n^„+i. Let us see how to do this. Again, we examplify the demonstration on a certain rank, 
but it is obviously valid for any. 
We have 

52 = ^Pi,j (55) 

i<j 

53 — ^ [pi,J and j+l,k+l ~ Pi, J and j+l,k] + Pi, J and j + l,fc,fc+l (56) 

i<j<k 

So, the passage from Sn to Sn+i is obtained by the formal linear operator: C — {. . ■)k+i — (• • ■)k + 
(■ ■ ■)k,k+i- Let us write now 5*2 in terms of C^^^ and C^"*: 

S2 = ^ [Pl and + l ~ Pi and + ^ W and i+lj j + l (57) 

i<j i<j 

The question now is: can we apply the operator C term by term to the equation H57() ? We have 

and i+l j + l] — Pi and i+1 j" + l,fc+l ^ Pi and i+l,j+l,k + Pl and i+l,j+l,k.k+l (58) 
^ [Pl and j+l j"] = Pi and i+l,j,k+l — Pi and i+l,j,k + Pi and i+l,j,k,k+l (59) 
[Pl and = Pi and i+l,j,j+l,k+l ~ Pi and i+l,j,j+l,k + Pi and i+l,j,j+l,k.k+l (60) 

On the other hand, the generic term of S3 (cf. eq H56|) ') can be expanded "inside" using the relations: 

Pi, J and j+l,k+l = Pi and i+l,j,j + l,k+l + Pi and i+lj,j+l,k+l (61) 

= Pi and — Pi and i+l,j,k+l + P/ and j+l (62) 

Pi, J and j"+l,fc = Pi and i+l,jj + l — P/ and i+l,j,k + P/ and i+l,j+l,k (63) 
P/,J and j"+l,fc,fc+l = P/ and i+l j" j + l,A;+l ^ P/ and i+l, j, j+l. k + Pi and i+l j + l,fe,fc+l 

— Pi and i+l j-,A;,fc+l + Pi and i+l,j ,j + l,/c,fc+l (64) 
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A rapid inspection shows that apphcation of C term by term in (|57|l gives the correct result. Moreover, 
one sees that = C(2) and £[(7(2)] = C^^) + C*^) and similar relations hold at any order: 



C 



in) 



N 



= C(n+i) + ^0;') 



(65) 
(66) 



A couple of polynomials (P„, Qn) can be associated to Sn- P-a represents the C terms and Qn the Cjv 
terms, with the mapping 

C(™) ^ x™. Thus, 



(Pi,Qi)-(X,0) 

{P2,Q2) = {X\X) 

{P,,Q,) = iX^ + X^,X^) 



{Pn+l, Qn+l) — (XPn + XQn, P„ 



(67) 
(68) 
(69) 
(70) 
(71) 



whence P„+i ~ XP„ + XP„^i. Symbolically, this recursion can be solved for P„, owing to Pi 

°2 



Pn 



X 



X + + 4:X\ X - ^/X^ + AX 



Qn is simply deduced from the formula (5„ = Pn-i- Symbolically, it yields 

Xt 



n=l 



l-X(r + r2) 



and 



n=l 



n=l 

From these results is derived a rewriting of G(/i) 

G(m) =1 + E(- + ^')"^^"^ + E(- + 



(72) 

(73) 
(74) 

(75) 
(76) 



n=l 



= 1 



T+1 



T+1 



(t + r^) ^ c,,j+i + (r + t2)2 ^ Ci^,+i jj+i + 
i •t<i 

(r + T^) ^ C,;,i+i,Ar + (t + r2)2 ^ Ci,i+ijj + i^ 



(77) 



The second sum can be slightly transformed : it can be shown (see Appendix) that 



(78) 



12 



whence 



C 



(n) \ ^ _ \ ^ 

N ~ / J '^jl,jl+l,j2,i2 + l,---,in,in + l / , CQ ,i 2 ,12 + 1 , • • • , + 1 , JV 

0^ii<i2< ■■<in 0<i2< --<4,i 



i3,...,i„,i„+l,iV 



(79) 



13<...<!„ 

(otherwise expHcitely stated, let us recall that the sums start at site 0). Whence 

,-2 1 \ ^2\2 



T + T 



71=1 



— i-^ ^(r + r2)«C(") + 1" + "^ 5:(r + r2)"Ci(„)^ 

(80) 



ri>0 



where 



As a result, one has 



Cl(n)JV — Ci_ij^^ij+i^i2,i2 + l,---,in,in + l. 

0<ii<i2<-<in 



(81) 



G(/i) = 



1+T 



1 + T + r2 



1 + ^(r + r2)"C(") + ^(r + t^)"^! 



(n)Af 



n^O 



(82) 



This rewriting is fruitful, since the sums are directly related to the Pfaffian theory: as explained in 
the Appendix, we can write 



g(m) 



1+T 



1 + r + t2 

2 



v/det[l + A{r],r]')C] with 



r/ = — T — r 



1 + r 



(83) 
(84) 

(85) 



(the definitions of A and C arc given in the Appendix). 
4.3 Asymptotic behaviour of G 

In the authors showed that the leading term of logdet[l + A{ri,if)C] remains unchanged if we put 
rj' = 0. In fact, in their case, a difficulty arises, due to the fact that the subdominant term (where 
rj' appears) diverges in a certain range of parameters and provides an analytic continuation of the 
persistence exponent. Here, the computation of the subdominant term cannot be performed exactly 
because the matrix AC does not have an exact continuous limit, but we think that it is precisely this 
fact that prevents the subdominant term to interfere in the analytic properties of the Idf. Anyway, we 
did not find any analytical problem in the Idf (see below) , and interpreted this as a "gentle" behaviour 
of the subdominant term. In the following the replacement ?/ = is implicitely made. 



Using logdct(A/) = Trlog(Af), we get 



logG(M) 



(-i)f 



Tr[(A(77,0)G)f] 



(86) 
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From the expressions of A and C, we have: 



Ti-[{A{r,,0)Cr] 



-(r + r^) 



E 



?1?2 ■ ■ -ip 



l^ii,...,ip^N 



(87) 



We get this first expression by assuming that the leading term is obtained when all indices are large 
(cf. [SI and see below). We assume temporarily that this term is O(logA^); thus we have necessarily 
for any a > 1 



E 



«1«2 ■ ■ - H 



loeN 



{ti + ilM + il) ■ ■ ■ i^l + ^l) log 



E 



lll2 .. .It, 



N^ii,...,ip^Na 



(88) 



log TV 
logo 



dui . . . dup 



(1 + e2("2-f i))(i + g2(u3-u2)^ . . . (1 + e2('"i""p)) 

(89) 



Let us now consider the large a limit in the integral and perform the change of variable vi = mi, 
V2 — U2 — Ml, . . .Vp = Up — Up-i (its Jacobian is 1). Defining s{x) — e^^/(l + e~^^) = l/2cosh(a;), we 
have 



logo 



dui . . . dup 



(1 + e2('"2-Mi)) . . . (1 + e2("i-"p)) 



oo /"log a 

dy I dvi I dv2 ■ ■ -dvp 6{v2 + - ■ ■+Vp~y)s{v2) . ■ ■ s{vp)i 



oo Jo 



log a X 



dk 
2^ 



dve'''''s{v) 



(90) 
(91) 



Before proceeding further, let us show that this term can be expressed differently, which will make the 
logA^ dependence explicit: 



E 



lll2 ■ ■ -Ip 



N 



ii,...A„ 



'E ^ E 



(h/ip) ■ ■ ■ (ip^i/ip) 



i{h/ipr + {^2/^pr)...{i^l/^pr + l) 



p\og{N) dxi . . . dxp^i—^ 
Jo {xf 



(92) 



X\ . . . Xp—\ 



(93) 

(because the general term of the ip sum is asymptotically equivalent to the integral times 1/ip)- This 
expression is however not as useful as the preceding, for it leads to an expression involving a solution 
of a Wicner-Hopf problem. 



Coming back to G(/i), eq. (|91ll and ^ yield 

logG(/i)~log(iV)x j —log 
'V2 



1 + -(t + t^) / dv e'^^'siv) 

TT ' 



5(m) 



-Arccos 



2^^ 



V2 



(94) 
(95) 



We verify that 5(0) = 0, -^'(0) = {E/\ogN) = 2/ti and,9"(0) = ^{-k-I) = {{E^) - {E)^) / log N , as it 
must be. The third derivative of g is 64(3 — 7r)/7r^ on zero. In the figure the aspect of /(C) is shown 
and its shape is compared with a parabola. The fluctuations above the mean value are comparatively 
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more probable than those below the mean value: intuitively a positive fluctuation of E regresses by 
annihilation pair by pair of the (excess) domain walls, and can be sustained as soon as the domain walls 
avoid each other in their erratic motion; on the contrary, a important negative fluctuation (i.e. E < (E)) 
will recede by energy refilling from the boundary, a process which is unavoidable if the population of 
domain walls is substantially lowered. This could explain qualitatively the counterclockwise tilt of the 
curve. 

As regards the right asymptotics of /(C), it is easily established (with the formula H94() which is the 
analytic continuation of g{fJ.)) that /(C) ~ — (ttC/S)^; this explains the vaguely parabolic look of /. 

Finally let us notice that the /(C) stops abruptly at C = 0; this means that the probability of 
observing an arbitrary small energy E decreases like N^^^^ (up to a 1 / \/\og N term); it illustrates the 
fact that an important "clearance" of the system due to an efficient (in dissipation) ffuctuation is not 
particularly unlikely. In |10| a same kind of Idf with a non diverging left branch was found in a very 
different context. 



-0.1 
-0.2 
-0.3 
S -0.4 
-0.5 
-0.6 
-0.7 



■ 0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 




Figure 4: The large deviation function /(C) for the Ising model (solid line) and the best parabolic 
approximation — (C — 2 / ir)'^ / {2g" (0)) (dots). The Idf does not continue in the range C < 0. 



4.4 Generalization to (an)isotropic Potts models 

The above calculations can be extended without difficulty to a g-states Potts model (driven by a "voter- 
model" dynamics: in [t, t + dt], each spin imitates its left/right neighbour with a probability dt and does 
nothing with probability 1 — 2dt) with an arbitrary input at ctq (the only requirements concerning this 
input is that it must be statistically stationary and have a finite time scale; we call this arbitrariness 
"anisotropic" for the different colours of the Potts model are not necessarily equiprobable). The energy 
is now defined as ej — 2(1 — (5ct ,o- +i) (we put the factor two to make the Ising model equivalent to the 
q = 2 Potts). If Po{ao) is the stationary probability of the spin (Tq, let us define 

b = 2x (^l-J2Poief^ (96) 

(6=1 in the Ising model; b — 2 corresponds to a A + A A model). 

Actually, this generalization does not modify deeply the preceding analysis : the reasoning on the 
coalescing walkers is unaffected by the multiplication of the "colours" , since either the walkers collapse, 
which implies automatically that they "carried" the same colour (remind that the walkers cannot be 
considered as domain walls or vice- versa), or they do not collapse, and in that case their colours are 
statistically determined by the ctq process, which can be considered as before as an ultrafast process. 
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It is easily shown that {E) ^ \og{N)2h/iT and similarly the function G is obtained from the preceding 
by making the replacement t br: 




Figure 5: Normalized Idf for different values of b for the anisotropic Potts models: 
/(C)/[&^/"((C))] 3.S a function of The Ising model corresponds to b = 1. 



The plots are 



the coordinates, so that their maximum is located at 2/tt and their curvature at the maximum is —1. 
Once normalized, these functions look very similar, and the physical informations are mainly contained 
in the first two moments: we have (£')/log7V - 2&/7r and {(E^) - (E)"^)/ \ogN ~ 4b[Tr + {tt - 4)b] / tt"^ . As 
a result, irrespective of the details of the injection, the anisotropic Potts models verify asymptotically 
the simple relation 



AE^ 



(7r-4) 



\oeN 



(98) 



A natural question would be to know if such a relation can still be verified if intrinsic dynamics, 
although dissipativc, is not equivalent to a voter model. 



5 Nonequilibrium "free energy" 

As it was possible to compute the large deviation function of E, we can ask whether a more precise 
description of the fluctuations of the energy field could be tractable. The first idea is to form uq groups 
of N consecutive spins and try to compute the probability to observe a certain coarse-grained energy 
profile {El, . . . ,Eng) when N is large. But this idea does not take into account the great intrinsic 
inhomogeneity of the dissipative systems: if the typical energy of the A'' first spins is O(logA^), the 
typical energy of the spins from indices A^ to 2N is only 0(1), which prevents to give a satisfactory 
description in terms of large deviation function. The second idea is to rectify by hand this inhomo- 
geneity, considering for instance variables like easy to show that they do not obey 
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a large deviation theorem in the large N limit. Actually the only way to construct a large deviation 
function capable of describing more precisely the energy fluctuations is to consider unequal blocks of 
spins: the first group gather spins ctq to a^v-i, the second spins from crjv to CTiv2_i,. . . , the last, spins 
from tr^no-i to aN"o-i- We are going to show that this inhomogeneous grouping of spins allows a 
relevant description in terms of Idf. 

For the sake of simplicity, we come back to the Ising model, but the generalization of the preceding 
paragraph could be processed along similar simple lines. As well, we forget for a moment the particular 
way to group the spins together and consider something more general: 

HifM) = (cxp J2 - <^^<^^+l)^ ^ (99) 

_ff is a characteristic function for M spins where fi is replaced by an arbitrary "field" /x. A multidi- 
mensional Laplace inversion performed over /x would give the probability of the energy density field. 

A similar procedure as for G leads to 

H{n) = l + ^npi + nTjPij + ■■■ (100) 

i i<j 

= 1 + Ti+T2 + ... (101) 

with Tj = (e~^^J — l)/2. The resummation procedure is more complicated than for G. Let us define 
Aj = Tj^i — Tj. The expansion of T2 using the recursion relation 1140|) leads to 

?2 = ^ TiAjPi and i+l.j + ^ nTjPi and i+ljj + l (102) 
i<j i<j 

Symbolically, the "r" term of Ti is transformed to a "tA + rr" term at the next step. A lengthy 
computation shows that 

^3 = ^ TiAj AkP I and i+l.j, k + ^ TjAjTkPl and i+l,j,k,k+l + 'y^.TiTkPl and i+l,k,k+l 
i<j<k i<j<k i<k 

+ ^ TiTj AkP I and i+l.j. j + l,k + ^ TiTjTkPl and i+l,j.j + l.k,k+l (103) 
i<j<k i<j<k 

=" tAA + tAt + TT2 + ttA + ttt" (104) 

where the last line is a symbolic rewriting of T^. Besides, the computation shows that (. . .)A is 
transformed into (. . .)(AA + At + T2)- Thus, as for G, these simple rules allow us to generate any 
order, which are rapidly increasing in complexity: 

T3 = f Af + f AA + fT2 + ff A + rff (105) 

^ T4 = tAtA + tAtt + rAAA + f AAt + f At2 + rf2 A + f 72? + rf At + rf AA + ttt2 + ttt A + tttt 

(106) 

These terms can be grouped into four categories: there are terms (called "u„" in the following) without 
any A, terms with a single and terminal A (called "vn"), terms like rAAr with multiple internals A 
but not ending with a A (called "?/;„"), and terms with multiple A ending with a A (called "zn"). 
Recursions relations exist between these fours sequences: one has for u and v: 

Ui = f , U2= TT, Un+1 = U„r + U„_lf2 (107) 

Vi = 0, Vn+1 = it„A (108) 

Let us recall that 
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• the elements t, ti and A of the symboUc algebra do not commute. 

• the element T2 is not equivalent to tt 
Formally, one can sum up the it„ and w„ terms: 

00 

U =^Un=T{\-T -T2T^ (109) 
n=\ 
oo 

V =^Vn=T{\-T-T2)-^'K (110) 
ri=l 

The w and z terms obey also recursion relations: 

Wn+l = VnT + Z„T + W„T + Z„^T2 (HI) 
Zn+1 ^ (Wn + Zn+Vn)A (112) 

where means "z„ without its final A". Defining W ~ X]n>3^n ^'^^ ^ ^ Sri>3-^"' ^ simple 
computation gives 

H{n) = l + U + V + W + Z (113) 

= 1+?(i-?-?2-A) ^ (114) 

This simple formal result hides a real complexity: let us recall that for instance 

rAfaA = ^ T,AjTlAec,^i+ij^k,k+i,e (115) 

i<j<k<e 

A useful rewriting of H is 

oo 

i/(/i) = 1 + [(1 - ? - fa)"' a] ''(l - ? - f2)-' (116) 

5.1 Recovering G 

If all the fij are constant, this sum is considerably simplified: in fact, all Aj but Am = t are zero, and 
T2 = TT such that 

H{n) = 1 + ?(1 - ? - ?2)"' + ?(1 - f - ?2)"' A (117) 
= 1 + r f;(r + r2)f + ^(r + r^)^ C^^+i' = G(/.) (118) 

G(/x) is recovered. 

5.2 Two groups of spins 

Let us now consider a situation with two clusters of fij: iij = j^i if j e [0, iV — 1] and fij = ^' if 
j £ [N, N"^ — 1]. In that case, two Aj only are nonzero {An = t — t' and A7V2 = r') and H is expressed 
as 

H{^i, ^i') = l + T{l-T- T2y^ +t{1-T- ?2)"'A(1 -T-T2y^ +t{1-T- ?2)"' A(1 - f - ?2)"'A 

(119) 
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We can write (/i, /i') as n') = H{0, i^i') + AH. It is easy to see that H{Q, fj,') (resp. AH) is made 
with all terms where the first of the smiis is r' (resp. t). Defining 

H{fl, ^i') = l + {T + T2){1 -T- T2r^ + {t + T2){1-T- f2)"^ A(1 -T- ?2)"^ 

+ {t + T2){1 -T- T2)-^A{1 -T- fa)"' A (120) 

we get 

H = 1 + {1 + t)AH + {1 + t')[H{0, fi') - 1] (121) 
^-t' + {1 + t)H + {t' -T)H{0,fi') (122) 

H is interesting for it is easier to compute than H: the last expression shows that once H is computed, 
so is H, since H{Q,fi') can be estimated along lines similar to those we followed for G. 

Let us define ^ = r + and ^' = t' + r'^. Similarly ^ (resp. ^') is the operator t + T2 applied on 
indices less (resp. greater) than N. As an example, 

^2^/ ^ ^2^/ ^ Ciiii+l,i2,i2+ljl (123) 

0^il<i2<N^ji 

Similarly Am and Ajy^ term the A operator where the index of application is precised. We have 

H{^,,|J') = l+ J2 E ?'^^?'+ E ?'?'^^^+ E e^N^'^N^ (124) 

Terms with a single A can be further simplified using (|79|1 : 

C^AatC''^ = C^'C'^Aat [ — C(p)_jv,(g-l),W2 + C(p),((}) - C(p-l),W,((j) + Cl,(p-2),Af,(g)] (125) 

(126) 

The definition of the symbols C... is generalized from (|8H) : C(p)_(g) deserves however a precision: it holds 
for a summation over ^ ii < . . . < < ^ ji < . . . < < iV^ (instead of . . . < A < ji < . . .). 
After a computation, we get 

i?(/i,/i') = ^^j^ E (e^r +^''Aw?'A^..) +^ Yl m^NC.^ipM) + ^N-^C,^ip),i,)M 

^ p+q>0 ^ P+q>0 

- Aat ^ f"' + r' ^ Aat^'^Ajvs + boundary terms (127) 

q^l q^l 

whence 

H(0, m') = i/(0, m')/(1 + ^') - E - + b-t. (128) 

and finally 

H{f,, ^^') = ]^ E (^''?' + ^A^^'Ajv^) + E ^''^"[^^C'i,(p),;v,(,) + Ajv.Ci,(p),(,),^.] + b.t. 

^ p+q>0 ^ P+q>0 

(129) 



19 



When T ^ t' , the formula for G is again recovered. It is interesting to note that H has here again a 
Pfaffian structure: if B is the skew-symmetric N'^xN'^ matrix where all elements but bi^j+i = — r, (ri + l), 
61, AT = — t^Aat, fei.Ar2 = — T^A^va, &Ar,Ar2 = — AArA^ra are zero, then 

1+T 



H{pi,^l') = - -v/det(l + BC) +b.t. (130) 

1 + T + T-^ 

By induction from the preceding case, we assume that the large deviation function can be computed 
as if Ajv = A7V2 = 0. In that case, we get 



(131) 
(132) 



where Bp is a trace term (we omit the p summation symbols) 

1 

2n 



Bl^^^l^ (133) 



5.3 General case 



For the general case, where no groups of spins are considered instead of two, the generalization is quite 
immediate and similar conclusions hold, that is, H is proportional to a Pfaffian (up to boundary terms) 
whose leading term is unchanged if off-second-diagonal terms are replaced by zero. 

As for the computation of G, there are two different ways to compute the leading term of Bp. For 
the inhomogeneous case, it is convenient to write 

B ^V^f V f £ {ii/ip){i2lip) ■■■{ip-i/ipyp 

'^^-'^ .^0 n.2,.^.-.=i ' ■ ' ([^lAd' + [*2Ap]2)([*2/zp]2 + [i,/ip?) . . . (1 + [^l/^p?)^l' 

(135) 

jV"o_i 



P X! ~ / "^^1 / dX2..- dXp-i £i^xi ■■■£.: 

i =0 *P "'0 "'0 -^0 



X\X2 ■ ■ ■ Xp-i 



'^^^-^ixi+xi)...{xi_, + m+xi) 

(136) 

Let us recall that we consider a particular profile where fii (i.e. is constant in each interval [0, N— 1], 
[N,N^-1], . . . , [7V"o-i,Ar»o_i], that is 6 = e(£^[log(i)/ log(7V)]). Thus, ^.,.1 = e(^[log(zp)/ log(7V) + 
log(a;i)/ log(iV)]) and £,i^xi — £,ip except for xi of order 1/N. As a result, in this case, 

N"-l ^P „1 „1 „1 XXX 

" ' £ ^ '0 ^0 ■ ■ ■ io '^-^ (.?+.i)';'(.^:r+"i)(i+^?) ^'''^ 

^p{iogNmor+£{'^r + ---+an-ir) t dx^ [\x2... f'dxp^i xix2...xp-i 



ixi+xi)...{xi^, + m+xi) 

(138) 



Thus, the coupling vanishes at the level of the large deviation function, and we have (the vector 
fi = (Ai(l)j ■ • ■ , m("-o)) holds here for the no different values of p for the no groups of spins) 



4=1 



Arccos 



V2 



(139) 
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Thus, it appears that the successive groups of spins are decorrelated at the level of the Idf. This 
decoupling shows that the situation is here completely different from those arising in conservative 
systems, and sheds light on the intimate relation between the conservative character of an observable 
and the appearance of severe correlations in a nonequilibrium stationary state. Here any fluctuation 
arising somewhere in the system affects only marginally its vicinity, for this fluctuation is mainly locally 
destroyed by the dissipation. 

6 Conclusion and Perspectives 

We exhibited an example of dissipative system in a stationary state where a global variable (the 
cumulative energy E) obeys a large deviation theorem; moreover its large deviation function f{x) 
appeared to be independent of the injection mechanism. This particular "boundary layer" structure 
made the exact computation of / tractable, and we were also able to compute the large deviation 
function associated with the probability of a particular energy profile. Due to strong inhomogeneitics in 
the stationary state, such a profile must be defined in a nonstandard way, and there is only one coherent 
way to do this. We proved also that the stationary state do not exhibit correlations as conservative 
systems do when pulled out of equilibrium, for the bulk dissipation prevents the fluctuations to be 
transported unchanged from one place to another. Incidentally, it is interesting to mention the reference 
|12| . where the authors come to quite different conclusions : they study the equilibrium properties 
of a system driven by Kawasaki+Glauber dynamics and show that long range correlations do exist, 
although the system be in equilibrium ; they conjecture thus that "correlations are more likely to be 
a generic feature of non reversible dynamics" . In our case however, where internal dynamics are also 
non reversible, correlations are not induced. Maybe the dynamical details play an important role, and 
nonlocal dynamics like Kawasaki's could impose correlations at the large deviation level. 

A natural extension of this work is to add a systematic drift to the system and test if it would 
be able to rebuild correlations ; besides, if we interpret it as a convection phenomenon in a realistic 
system, the study of such a competition (and balance) between the dissipation and the transport 
in the structuration of the fluctuations could give interesting insights toward the characterization of 
dissipative NESS. 

But some care must be taken in doing a modification of internal dynamics of a dissipative system: 
actually, it is quite easy to destroy the large deviation "structure" of the energy profile. In fact, the 
appearance of strong correlations in a stationary state seems to be slightly antithetical with a large 
deviation theorem, since the latter is mathematically expressed for a sum of N independent random 
variables; if we know that as a rule weak correlations do not prevent a large deviation theorem to 
hold, it is obvious that strongly correlated systems have no reason to verify this theorem anymore. In 
this respect, the results of ^ seem at first sight paradoxical, since strong correlations are present in 
their systems; but actually, if their correlations are strong in the sense they "connect" arbitrary distant 
places in the system, they are nevertheless weak as regards their intensity, since they are inversely 
proportional to the number of particles. 

7 Appendix 

The central point of this model is that probabilities c..., whatever the event symbolised by the . . ., can 
be expressed in terms of the only. To see this, let us consider the recursion relation verifies: 




(140) 

(141) 
(142) 



coj = 1 V j ^ 
CiA ~ V i ^ 1 
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The equation (|140|1 comes from an explicitation of the next move of the walkers: there are 4 equiprobable 
moves and for instance Prob[(z, j)|next move is (i ^ i + 1)] = Q+ij. The boundary conditions 1] 142(1 
ensure that this relation is valid even if i + 1 = j. Similarly the b.c. ((141|l plays the same role for the 
event i i — 1 when i = I. 

Very similar relations can be written for any Cij^k,... (the r.h.s of the recursion relation has 6 terms 
for Cij_k, 8 terms for Ci_j^k,e, etc. . . ). It is then just a matter of calculation to verify that the solution 
of the recursion relation for Ci_j^k and Cij^k,i are 



These relations can be generalized to higher orders: 



-ll,!2,---,l2n + l 
Cii,i2,---,'j2n 



^ll,...,t2n 



T > e(a)cj j Cj ■ i 

' o-eS2„ 



-"'2,---,*2re + l 



^2n — l "^271 



(143) 
(144) 



(145) 
(146) 



In the last expression. S2n is the set of permutations of {1, 2, ... , 2n} and the Cij are antisymmetrised: 



-Ci,j if j > i- The formula (|146|l show that the Qj. 



are Pfaffians associated with the 



antisymmetric matrix C = (cij): as demonstrated in if A is the N x N matrix 



/ 1 c \ 
-10 10 

0-101 
. 

. 1 

V -e -10/ 



(147) 



We have the relation 



det (Id + AC) 



.N 



-S. yici.N ^v^Yl ci,is+i. 



(148) 
(149) 
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